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White dwarfs are almost completely degenerate objects that cannot obtain energy from the 
thermonuclear sources and their evolution is just a gravothermal process of cooling. The 
simplicity of these objects, the fact that the physical inputs necessary to understand them 
are well identified, although not always well understood, and the impressive observational 
background about white dwarfs make them the most well studied Galactic population. 
These characteristics allow to use them as laboratories to test new ideas of physics. In this 
contribution we discuss the robustness of the method and its application to the axion case. 

1 Introduction 

White dwarfs are almost completely degenerate objects that cannot obtain energy from the 
thermonuclear sources and their evolution is just a gravothermal process of cooling. Globally, 
the evolution of the luminosity of a white dwarf can be written as: 

/■Mwd j T j-M wn /dP\ dV 

L + L„ + L x = - c v —dm- T(—) —dm+(l s + e s )M s + e x (1) 

J dt J Q \dTJv,X dt 

where the l.h.s. of this equation represents the sinks of energy, photons, neutrinos and any 
additional exotic term, while the r.h.s. contains the sources of energy, the heat capacity of the 
star, the work due to the change of volume, the contribution of the latent heat and gravitationl 
settling upon crystallization (the term M s is the rate of crystallization) and, finally, the last 
term represents any additional exotic source of energy pQ . There are two ways to test the theory 
of white dwarf evolution, namely, studying the white dwarf luminosity function and using the 
secular drift of the pulsation period of variable white dwarfs. 

During the cooling process, white dwarfs experience some phases of pulsational instability 
powered by the k- and the convective driven- mechanisms [2]. Depending on the composition 
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of the atmosphere variable white dwarfs are known as DOV, DBV and DAV. These stars also 
known as PG1159 or GW Vir stars, V777 He stars and ZZ Ceti stars respectively. Variable 
white dwarfs of different compositions occupy different regions in the Hertzsprung-Russcll dia- 
gram. The value of the pulsation period indicates that these objects are experiencing g— mode 
non-radial pulsations, where the main restoring force is gravity. One of the main characteris- 
tics of these pulsations is that they experience a secular drift caused by the evolution of the 
temperature and radius. For qualitative purposes this drift can be approximated by [3]: 

dlnn dlnT dhxR , , 

— ^— ^-a—— + b— — 2 
dt dt dt w 

where a and b are constants of the order of unity that depend on the details of the model, 
and R and T are the stellar radius and the temperature at the region of period formation, 
respectively. This equation reflects the fact that, as the star cools down, the degeneracy of 
the plasma increases, the buoyancy decreases, the Brunt- Valsala frequency becomes smaller 
and, as a consequence, the spectrum of pulsations gradually shifts to lower frequencies. At the 
same time, since the star contracts, the radius decreases and the frequency tends to increase. 
In general, DAV and DBV stars are already so cool (and degenerate) that the radial term is 
negligible and the change of the period of pulsation can be directly related to the change in 
the core temperature of the star. The timescales involved are of the order of ~ lfV 11 s/s for 
DOVs, - 1CT 13 to - 1CT 14 s/s for DBVs and - 1CT 15 to - 1(T 16 s/s for DAVs. Therefore, 
the measurement of such drifts provides an effective method to test the theory of cooling white 
dwarfs. This measurement is a difficult but feasible task, as it has been already proved in the 
case of G117-B15A, a ZZ Ceti star [4]. These properties allow to build a simple relationship [5j[6] 
to estimate the influence of an extra sink of energy, axions for instance, on the period drift of 
variable white dwarfs: (Xx/^modei) ~ (n o b s /n mo dci) — 1 where the suffix "model" refers to 
those models built using standard physics. 

The white dwarf luminosity function is defined as the number density of white dwarfs of a 
given luminosity per unit of magnitude interval: 

n(l)= $(M)*(t)T cool (l,M)dM (3) 

J Mi 

where t satisfies the condition t = T — t coo \(l,M) — tps(M) and / = — log(L/L Q ), M is the 
mass of the parent star (for convenience all white dwarfs are labeled with the mass of the 
main sequence progenitor), i coo i is the cooling time down to luminosity I, t coo \ = dt/dM\, \ is 
the characteristic cooling time, M s and Mi are the maximum and the minimum masses of the 
main sequence stars able to produce a white dwarf of luminosity I, tps is the lifetime of the 
progenitor of the white dwarf, and T is the age of the population under study. The remaining 
quantities, the initial mass function, <£(M), and the star formation rate, ^(t), are not known 
a priori and depend on the properties of the stellar population under study. The computed 
luminosity function is usually normalized to the bin with the smallest error bar, traditionally 
the one with I = 3, in order to compare theory with observations. Equation ([3]) shows that in 
order to use the luminosity function as a physical laboratory it is necessary to have good enough 
observational data and to know the galactic properties that are used in this equation (the star 
formation rate, the initial mass function and the age of the Galaxy). Fortunately, the bright 
branch of the luminosity function is only sensitive to the average characteristic cooling time of 
white dwarfs at the corresponding luminosity when the function is properly normalized. The 
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reason is twofold [7], on one hand the stellar population is dominated by low mass stars and 
on the other the lifetime of stars increases very sharply when the mass decreases. The result 
is that the old galactic populations arc still producing bright white dwarfs and the number of 
such stars at each luminosity bin is the sum of the contributions of all the different episodes of 
star formation. 
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Figure 1: Power emitted in the form of photons (continuous line), neutrinos (dotted line) and 
axions (dashed line, arbitrary <7 a cc) by a 0.6 M Q white dwarf versus its bolometric magnitude, 
a function that monotonically increases with time. The cross represents the position of G115- 
B15A. 



2 The axion case 

We note that in white dwarf case only DFSZ axions have to be taken into account, since 
electron bremsstrahlung is the dominant process in white dwarf interiors. The axion emission 
rate under these conditions is given by [5]: e ax = 1.08 x 10 23 (gl cc /4ir)(Z 2 A) T 7 4 F(T, p) erg/g/s 
where F takes into account the Coulomb plasma effects, Tj is the temperature in units of 10 7 
K, Z and A are the atomic and mass numbers of the plasma components, respectively, and g aoc 
is the strength of the axion-electron Yukawa coupling. Figure [T] shows the energy losses of a 
typical white dwarf star. Since the neutrino emission is e„ cx T s and L p hot oc T 2 7 the luminosity 
function allows to disentangle the contribution of the different mechanisms of energy loss. When 
this method is applied to the luminosity function of white dwarfs the value of <? aee that best fits 
the luminosity function is 1.1 x 10~ 13 , but variations of a factor two are still compatible with 
the observations [9J, [7] . Finally, the most recent analysis of G117-B15A shows that the value 
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of g a ee quoted here is compatible with the secular drift of the pulsation period, which gives 
support to the necessity to include an extra cooling term in the white dwarf models 10 . 

3 Conclusions 

There are two independent evidences, the luminosity function and the secular drift of DAV 
white dwarfs, that these stars are cooling down more rapidly than predicted. The introduction 
of an additional sink of energy linked to the interaction of electrons with a light boson (axion, 
ALP,...) with an strength g acc ~ 10~ 13 solves the problem satisfactorily. Of course, the 
remaining uncertainties, both observational and theoretical, still prevent to claim the existence 
of such interaction. A direct detection under laboratory conditions, like those of the CAST 
helioscope QTj or the shining through the wall experiments |12j . could provide unambiguous 
evidences. In this sense, if the current result is interpreted as due to axions, such particles 
should have a mass, m a ~ meV and should be coupled with photons with g a7 ~ 10 -12 GeV -1 . 
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